Convergence of a renormalization group approach to dimer-dimer scattering 
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We study the convergence of a functional renormalisation group technique by looking at the ratio 
between the fermion-fermion scattering length and the dimer-dimer scattering length for a system 
of nonrelativistic fermions. We find that in a systematic expansion in powers of the fields there is a 
rapid convergence of the result that agrees with know exact results. 
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I. INTRODUCTION 

The atomic physics of ultra-cold Fermi gases is one of 
the places where we can make a detailed link few-particle 
and many-body physics. In fermionic systems we can 
trace the effect of an attractive force from the bound 
states or resonances in few-body systems to the pairing 
occurring in the many-body sector. At low energy and in 
cold gases this physics is governed by a single parameter, 
the central S-wave scattering length ap, determining the 
scattering at threshold. We can also tune the value of 
this scattering length making use of recent advances us- 
ing Feshbach resonances. For negative scattering length 
the gas is in the weak-coupling BCS state. For positive 
values of op bound states of two fermions- "dimers" -form 
and these can lead to a Bose-Einstein condensate (BEC) 
The size of dimers is determined by the fermion- 
fermion scattering length and their binding energy is of 
order l/ap. 

If we concentrate on the case of deeply-bound dimers, 
then for a sufficiently dilute and cold gas of dimers the 
main dynamical quantity characterising their interaction 
is now the dimer-dimer scattering length ob, which is 
a induced by the fermionic scattering. The exact rela- 
tion between dimer-dimer and fermion-fermion scattering 
lengths a,B = 0.6of was established in Ref. Q by solving 
the Schrdinger equation for two composite bosons inter- 
acting with an attractive zero-range potential. Unfortu- 
nately, this method is difficult to extend to the many- 
body case. Therefore, it is useful to study the ratio 
clb/cif m an approach which can be used both for few 
and many-body problems. 

One such method is the functional or "exact" renor- 
malisation group (ERG) for the one-particle irreducible 
effective action, which is the approach applied in this 
paper to calculate as- (For reviews, see Refs. 0, |U). 
This technique has been previously used to study a va- 
riety of physical systems, from systems of nonrelativistic 
fermions [B4H| and bosons [l2j to quark models [l3| and 
gauge theories [13] ■ It is based on the scale-dependent 
quantum effective action Tk, where the action at scale k 
contains the effects of field fluctuations with momenta q 
larger than k only. In the limit k — > + all fluctuations 



are included and the full effective action is recovered. In 
practice the scale dependence is introduced by a set of 
fc-dependent cutoff functions R(q), which suppress the 
effect of modes with q < k in the path integral for the 
action by giving them a large k dependent mass. The 
minimal conditions satisfied by the functions R(q) are 
then that they should vanish in the limit k — > + and 
scale like k a with a > when k — > oo . 

With this prescription the average effective action at 
very large k converges to the classical action of the 
theory — here nonrelativistic fermions with a local zero- 
range interaction. Since one can show that we get the 
effective action for the theory as k — > + , the endpoint at 
k = of the exact solution of the functional RG equation 
should be independent of the choice of cutoff. However, 
in practice, truncations of the action inevitably lead to 
some cutoff dependence of the results. We can use this 
dependence as one possible measure of the quality of the 
truncation. With this tool, we shall analyse the conver- 
gence of a low-energy long-wavelength expansion as we 
increase the complexity of the many-body truncation in 
the effective action. 

As discussed in great detail in the literature, e.g. 
Ref. [HI, Chapter 16], the idea of the (quantum) effec- 
tive action T is to introduce an object that generalises 
the concept of the classical action of a field theory to 
include all quantum effects, but still depends on a set 
of classical external fields (the dual of the usual exter- 
nal sources in the partition function). The ground state 
is the minimum of the action, and we can generate all 
Green functions by appropriate derivatives of T. Apart 
from for very simple models, it is very hard to explic- 
itly evaluate the full effective action, which is given by a 
complicated path integral. The ERG technique used here 
gives a scale-dependent mass-gap to the low-momentum 
modes in the path integral, so that their fluctuations are 
suppressed. The action then flows from the trivial case 
of all fluctuations suppressed (scale k is infinite) to the 
full quantum effective action for zero scale. Since this 
approach in principle incorporates the full quantum dy- 
namics of the underlying theory, it is called "exact" . 

The flow of the effective action satisfies the functional 
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differential equation Q 
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where I^ 2 ' is the second functional derivative with re- 
spect to the fields, and the cutoff functions in the mass- 
like term R(k) drive the RG evolution. The opera- 
tion STr denotes the supertrace fl6| taken over energy- 
momentum variables and internal indices and is defined 
by 
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The evolution equation for the average effective ac- 
tion thus has a one-loop structure, but contains a fully 
dressed, scale-dependent propagator (r^ 2 ^ — Thus, 
despite its apparently simple form, Eq. (TJJ) is actually a 
complex functional differential equation. In the absence 
of general methods to numerically solve such equations 
we must resort to approximations. One common ap- 
proach is to parametrise the effective action with a small 
number of terms, which turns the evolution equation into 
a system of coupled ordinary differential equations for the 
numerical coefficient of each term. These equations can 
then be solved numerically. Here we study possible trun- 
cations for fcrmionic few-body systems, and our choice of 
ansatz for the action is motivated by both ERG studies 
of many-body systems 0] and effective field theories 
(EFTs) for few-fermion systems [l7|. The technique we 
use is similar to the one used in Ref. @ and Ref. @. 
The parametrisation of the effective action is given by a 
local expansion in fundamental and induced fields, with 
a simple gradient expansion for non-local terms. 

There are a number of questions about the approxi- 
mations made in writing down the action. We shall try 
and answer two of these in this paper: Firstly, does the 
low-energy truncation, where we expand to first order in 
energy and to first order in field-gradient squared con- 
verge? Secondly, are these results independent of the 
choice of cut-off function? These two questions are not 
independent: If we have a converged truncation of the ef- 
fective action, we must have independence of the cut-off 
function. Reversing this, we can hope that the depen- 
dence on the cut-off can be used as one of the signals of 
convergence. 

The renormalisation functions R(q) is not completely 
arbitrary: as the cutoff scale tends to infinity, we de- 
mand that the action be purely fermionic, containing a 
contact two-body interaction without derivatives. Such 
interactions are of interest also for interacting nucleons 
In both these cases, it is convenient to add an aux- 
iliary composite boson field (which is really the dimer 
field) to the problem by making a Hubbard-Stratonovich 
transformation. This then replaces the two-body inter- 
action by a Yukawa-type coupling between the fermions 
and the auxiliary boson. At the extreme end of the scale, 
the boson is not dynamic (there is no energy or momen- 



tum dependence of this field in the action), but a kinetic 
term for the boson is generated in the RG evolution. 



II. BUILDING BLOCKS 

We define the low-energy effective action following pre- 
vious work (see (B^SKIH) as an expansion in the fermion 
fields i/j—a, two-component spinor field, where the index 
denotes fermion spin- and the scalar dimer boson field 
4> obtained from bosonising the bilinear fermion field 
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In contrast with the many-body case, in the few-body 
sector there is no wavefunction and mass renormalisa- 
tion of the fermion field. The quantity Tl(x,x', k) is the 
scale-dependent boson self-energy and U2 parametrises 
the boson-boson interaction which is generated by the 
evolution. The latter is equivalent to a four-body inter- 
action in terms of the underlying fermions. The term 
proportional to A describes the fermion-dimer scattering 
(three body in the fermions). The final two terms are 
four-body again: the term proportional to g' describes 
the scattering of two dimers, where one dimer breaks 
up into two fermions, and finally the v term describes 
the scattering of two fermions from a dimer, without 
changing character. We have not included the local two- 
fermion interaction, which has zero coupling constant due 
to our choice of starting model at large (infinite scale), 
where only g is nonzero. 

The evolution of the boson self-energy is given by 



d k n(x,x',k) = 



6<f>(x')S^(x) 
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although from now on we shall express all evolution in 
momentum space. Note that only fermion loops con- 
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tribute to the evolution of the boson self-energy in vac- 
uum. 

Using a gradient expansion of the action, we can define 
boson wave- function and mass renormalisation factors by 
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where £p> = —l/(Ma F ) denotes the bound-state energy 
of a pair of fermions. It is relatively straightforward to 
show that Z m (0) — Z<f,(0), independent of the cut-off 
function. There are very few cut-off functions that pre- 
serve Galilean invariance to quadratic order in momenta, 
but one favourite sharp cut-off 19] 



RF(q,k) = ^fe(k-q), 



(7) 



can be shown to preserve the identity Z<j,(k) = Z m (k). 
Such sharp cut-off functions are difficult to apply in 
medium, since the evolution of Z m will contain ambigu- 
ous terms containing S functions and their derivatives 
arising from the first and second derivative of the cut-off 
function. Such difficulties can be bypassed here since II 
can be evaluated directly. We first impose the boundary 
condition that the scattering amplitude in the physical 
limit k — > reproduces the fermion-fermion scattering 
length, 
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Here Po (P) denotes the total energy (momentum) flow- 
ing through the system and p = y / 2MPa — P 2 /2 is the 
relative momentum of the two fermions. Integrating the 
resulting ERG equation gives @ 
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This shows that Galilean invariance is preserved only to 
lowest order. 

The complexity of sharp cut-off functions in medium 
implies that it is of interest to study smooth cut-offs, 
where both Z's must now be calculated independently. 
A suitable set of smooth functions can be parametrised 
as 27] 
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A. Mean-field 

In vacuum, it is easy to show that 
^ r d 3 q 1 
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The evolution of the boson-boson scattering amplitude 
follows from 



-d k u 2 {£ Dl k) = s . x2lc _ n ^z_ Li2rc _ m d k T\,p = o. 

(11) 

The driving term of this equation can be separated into 
fermionic and bosonic contributions containing dkR F and 
dhRs, respectively. We first look at the "mean-field" 
result, where bosonic contributions are neglected. The 
evolution of u 2 is then given by 



dkU 2 
where 
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Equation (|12p is integrable, and we find with w 2 (oo) = 
that 



u 2 (k) 
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u 2 (0) = —M 3 g 4 a F . 
The scattering amplitude at threshold is 



(14) 
(15) 
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giving the well-known mean- field result a R = 2a F [25| 
which is far from the exact value as = 0-6a F 0. This im- 
plies that beyond-mean-field effects such as dimer-dimer 
rescattering must be considered. 



B. Boson and mixed loops 

To include such effects we must take into account the 
loops containing boson propagators. After some algebra, 
we find the driving term 
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where 

E BR ( q ,k) = ^l q 2 +u 1 (k)+R B (q,k) (18) 

and 

u x {k) = -U{£ D ,0,k). (19) 

In the case of a sharp cutoff, we choose the bosonic cutoff 
function to be as close as possible in form to the fermionic 
one, 

R B (q, k) = Z m (k) {CBk ^ M ^ 0{c B k - q), (20) 

apart from the addition of a parameter c B , which sets 
the relative scale of the fermionic and bosonic regula- 
tors, and a factor of Z m . The latter has the important 
advantage of leading to a consistent scaling behaviour, 
so that all contributions to a single evolution equation 
decay with the same power of k for large k. Moreover 
it also gives a^-scaling, where all terms in a single equa- 
tion have the same dependence onap. Neither of these 
two conditions is actually required, and we shall show 
below that removing such restrictions may give us access 
to interesting information. 

It has been shown for bosonic systems [l2j , and it also 
shown in detail in Fig.[2]below, that the simple two-body 
truncation is insufficient. The complete set of contact in- 
teractions in the four-body sector includes not just the 



dimer-dimer term it2, but also terms where the dimer is 
broken into two fermions, Eq. ([3]). The evolution equa- 
tion for any of the couplings is described by an expansion 
about the energy of the bound state pole for bosons, and 
half that energy for fermions, e.g., 

S 4 STv[d k R(T^ -R)- 1 } 

k 2 6<f>\£ D ,0)8<t>(£ D ,0)S^(£ D /2,0)5i;(£ D /2,0y 

(21) 

The technique to evaluate such contributions can most 
compactly be written as a combination of a diagrammatic 
and algebraic approach. We first evaluate the skeleton 
diagrams that contribute to a given vertex. There are 
three distinct contributions to the running of A, com- 
ing from ladder, triangle and box diagrams, as shown in 
Fig. [T] For the other couplings, we have a large num- 
ber of diagrams that can contribute. The most difficult 
calculation is for g', where we have loops with up to six 
internal lines. Each diagram involves a single integration 
over an internal four momentum. Since our cut-off func- 
tion is only momentum-dependent, i..e., is independent 
of the energy variable, we can perform the energy integra- 
tion by a contour integration, enumerating the poles by 
solving linear algebraic equations. The insertion of the 
k derivative of the cut-off function on each leg can then 
be achieved afterwards, by a functional derivative of the 
resulting integrals with respect to the cut-off function. 

The resulting equations can be written in a compact 
form as, 
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III. RESULTS 



ically independent of the starting scale provided this is 



A. Full evolution 



In theory we should start integration of the evolution 
equations at infinity-in practice the results are numer- 
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FIG. 1: [Colour online] The skeletons of the diagrams that 
contribute to the evolution of the couplings. Dashed lines 
represent bosons, solid lines fermions. The dots represent 
the interaction vertices. In a) we show the diagrams for the 
evolution of ui, Z m and Z$ in b) for U2, c) for A, d) for g' 
and e) for v. 



chosen to be at least hap ~ 100. For k 1/ap the 
system is in the universal "scaling regime" , and a sta- 
ble fixed point, see below, governs the evolution until k 
becomes comparable with 1/ap. 

In this scaling regime, we can determine the behaviour 
most easily in terms of dimensionless "scaling variable" 
k = kap by defining the four dimensionless functions 

Cj(«) 

3„3 -3, 



u 2 (k) -4 g Nrapn'^coin), 
A(k) -4 g 2 Ma F n~ 2 c 1 (K), 

g'{n) -4 g 3 M 2 a%K- 4 c 2 (n), 
v(k) -4 g 2 Ma|^~ 5 c 3 (K). 



(23) 



The functions q satisfy a set of dimensionless differential 
equations. For each level of truncation of the effective 
action we can determine the resulting nontrivial fixed 
points, and the evolution close to these points, as given 
by the anomalous dimensions r\j , 



- Jp 
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For each truncation considered here we find only one sta- 
ble fixed point, see Table fl] There we give the value for 





FIG. 2: [Colour online] Ratio of dimer-dimer to fermion- 
fermion scattering lengths as a function of the relative scale 
parameter cb- The black (dash-dotted) curve shows results 
for the minimal action; the red (dashed) curve shows the ef- 
fect of adding the local three-body term; The green (solid) 
curve displays the full local four-body calculation, and is en- 
larged in the inset. The purple (dotted) line shows the effect 
of using a smooth cut-off (a — 0.5) on the full results. 



fixed point and anomalous dimensions for the sharp cut- 
off (O with the parameter eg = 1 in the bosonic cut-off 
([20| . These results are somewhat dependent on eg, as 
expected, and are also not totally cut-off independent, 
but seem to be reasonably stable under perturbations. 
The dimension for A, 
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IT 



3.10355 



(24) 



should be compared to the exact result 4.32244, found by 
Griesshammer and others [20-22]. The lowest anomalous 
dimension for the four-fermion sector, 4.19149, should be 
compared with the value 5.0184 obtained numerically by 
Stecher and Greene [23[ (see also Ref. (24|). 

If we start the evolution from the stable fixed point, 
we can then carefully trace this back to finite k. The 
behaviour of as jap as a function of cb for both a sharp 
and a smooth cut-off with a = 0.5, Eqs. and (|10p . 
is presented in Fig. [51 As we increase the complexity of 
the truncation, we see a rapid convergence: the inclusion 
of A reduces the size of as by a factor of about 2 (for 
small cb), and adding the remaining terms in the action 
pushes the result down even further to agree with the 
Schrdinger equation results; in that case we also see a 
very weak dependence ones in the region around 1. 

Note that in all cases the dominant contributions to 
o-b/clf for large cb come from the boson- loop terms in 
the equation for u%. Since these do not depend on the 
three-body coupling A, the curves approach each other. 
Moreover, this limit corresponds to integrating out the 
fermions first, which generates a non-zero value for u 2 at 
the start of the bosonic integration. In the limit cb — > oo, 
this coupling is driven to the trivial fixed point of the 
RG equations, u 2 = 0, since we have no terms to cancel 
the linearly divergent boson-boson loop diagram and the 
diagrams with three-body and four-body couplings are 
all too weak to alter this behaviour. 
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TABLE I: The coefficients and the eigenvalues of the stable fixed point for various level of truncation 
truncation fixed point anomalous dimensions 



u 2 c = 0.0656565 

u 2 , A co = 0.037878, ci = -0.322441 

U2,X,g',v c = 0.0356522, Ci = -0.322441, 
c 2 = 0.102042, c 3 = -10.488 



3.1728 
3.09969, 3.10355 
6.10943, 5.15046, 4.19149, 3.10355 



On the other hand, the main contributions for small 
cb come not only from the fermion loops, but also 
from mixed fermion-boson loops, which appear in the 
equations for the many-body couplings. In particular, 
the mixed boson-fermion loop diagrams containing the 
fcrmionic cut-off contribute to the evolution of A, v and 
g' , even when the bosonic degrees of freedom have been 
integrated out. As a result, inclusion of the three-body 
term A already leads to a significant deviation from the 
mean-field result, clb/clf = 2, that persists in the limit 
cb — > 0. With g' and v included we see convergence 
close to the exact result, even when cb = 0. We seem 
to have approximate convergence for a range of values of 
cb, probably best near eg = 1, but we can probably use 
any cb < 1.5. The strange results obtained for very large 
values of cb should not be taken too seriously, since they 
are based on an incorrect approach: we make the induced 
bosonic degrees of freedom dominate in the early stages 
of evolution (large k) . The other extreme, albeit naively 
equally incorrect, actually seems to produce sensible re- 
sults. This corresponds to freezing the renormalisation 
of the bosonic degrees of freedom, while still allowing an 
evolution of the coupling constants driven by the evolu- 
tion of the fundamental fermionic fields, suggesting that 
this may be a sensible and simplifying approach. This 
may have important practical consequences for calculat- 
ing in the many-body system: If we can integrate out the 
bosons at every stage, and only let the fermions evolve, 
the calculations become much simpler. 

Arguments based on "optimisation" of the cut-off func- 
tion, see Ref. (26[, indicate that one should choose the 
cut-off to try to maximise the rate of convergence for our 
expansion of the action. In this case there is a stationary 
point for cb close to 1, which agrees with the natural 
assumption that cb = 1, where bosons and fermions are 
renormalised at the same rate is the optimal choice. 

To provide a further check of convergence, it pays to 
use a different form for the boson regulator, and pur- 
posefully violate both the uniform a^-dependence and 
the uniform scaling for large k. If we have convergence, 
the results should remain independent of the cut-off, and 
thus also independent of ap and the shape of the cut- 
off function.. The simplest form we can choose uses the 
smooth cut-off function (fTU)) . where we renormalise the 
bosons as (note the absence of Z m ) 

Rb = ^ Cg ^ 9o(q,kc B )- 

IVtlB 

With this choice we expect the ratio qb/cif to be a 




FIG. 3: [Colour online] Ratio of dimer-dimer to fermion- 
fermion scattering lengths as a function of the relative scale 
parameter cb. The black (right-slanted hash) curve shows 
results for the minimal action; the red (vertical hash) curve 
shows the effect of adding the local three-body term; The 
green (left-leaning hash) curve displays the full local four- 
body calculation. The width of each curve denotes the sensi- 
tivity of the result to changes in at af = gM/h 2 



function of ap] we have first performed calculations at 
Op = 1, and looked at the sensitivity to changes in aj?. In 
Fig. [3] we see a strong dependence onap for any value of 
cb^O for the two-body truncation, a weaker dependence 
for the three-body case, and a very weak dependence for 
the four-body long as we consider cb < 1. For 

large cb > 2, where we let the bosons dominate, we have 
universally poor result whatever the truncation, in agree- 
ment with the discussion above. 

The convergence for the full four-body truncation can 
be seen even more clearly in Fig. |4l where we show a 
contour plot of the ratio of dimer-dimer and fermion- 
fermion scattering lengths as a function of oj? and cb- 
We note the excellent convergence in a large region of 
the parameter space, with allows a conservative estimate 
o,b/clf — 0.58 ± 0.02, in excellent agreement with the 
exact result. 



IV. DISCUSSION 

We have gathered considerable evidence for the con- 
vergence of a gradient expansion of the quantum effec- 
tive action for a system of a few dilute fermions inter- 
acting through a pairwise attractive force. The resulting 
dimers, fermion-fermion bound states, scatter in the way 
predicted by exact calculations, if we expand the polar- 
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FIG. 4: [Colour online] Ratio of dimer-dimer to fermion- 
fermion scattering lengths as a function of the relative scale 
parameter cb, and the logarithm of the scattering length of 
(divided by gM/h 2 ). Notice the insensitivity to clf for cb < 1, 
a clear indication of convergence. 



glected non-locality of all but the simplest terms, and 
have only included time derivatives to match the mo- 
mentum dependence. In principle, we can add any 
momentum-dependence to any terms without problems; 
it is much more difficult to add energy dependent terms; 
with a momentum-dependent cut-off we are limited to 
first order terms only. Fortunately, it appears that we do 
not need such complications! As long as we study low- 
energy physics, which is exactly the situation here and in 
the many-body situation, that is not too surprising. 

What does come as a surprise is that we seem to be 
able to fix the bosonic fields, and have an RG flow driven 
by the evolution of the fermionic fields only, while still 
obtaining good results. This is probably due to the fact 
that the evolution of the induced dimer degrees of free- 
dom can be thought of as driven from the basic fermionic 
degrees of freedom through the coupling constants. This 
requires confirmation for finite density many-body sys- 
tems. 



isation to second order in momenta, and include all the 
local four-body terms in the fermion and dimer fields. 

Of course the expansion is only complete in terms of 
the number of fields that enter the action: we have ne- 
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